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SOLUTION  OF  MAXWELL’S  EQUATIONS  IN  A  MAGNETO-IONIC  MEDIUM  WITH  SOURCES 


1,  Introduction 

The  knowledge  of  field  solutions  of  Maxwell’s  Equations  with  source  terms 

is  of  considerable  importance  for  many  applications.  In  this  paper  we  have 

derived  a  complete  solution  of  the  field  equations  in  a  form  which  is  suitable 

for  the  calculation  of  the  fields  of  arbitrary  distances  from  the  source. 

For  an  infinitesimal  dipole,  the  very  near  fields  are  found  in  a  closed  form 

involving  trigonometric  functions.  The  remainder  of  the  solution  is  expressed 

in  terms  of  a  finite  range  integral^  the  integrand  of  which  is  finite  everywhere. 

Hence,  the  solution  is  very  suitable  for  numerical  calculations.  It  is  also 

pertinent  to  point  out  that  the  solutions  derived  hold  for  a  very  general  nature 

of  the  elements  of  the  €  tensor;  for  instance,  the  solutions  are  valid  for 

complex  elements  and  hence  are  useful  for  the  treatment  of  a  medium  with  finite 

losses.  Two  other  previous  contributions  on  this  topic,  e.g.  by  Bunkin1  and 
2 

Kogelnik  do  not  give  a  general  solution  such  as  the  one  derived  in  this  paper. 

2.  Derivation  of  Matrix  Equations  in  a  Magneto-Ionic  Medium  with  Sources 

In  this  section  we  shall  derive  the  matrix  equations  for  the  fields  in  a 

homogeneous  medium  with  tensor  dielectric  properties,  for  the  case  of  impressed 

electric  current  sources.  The  equations  will  subsequently  be  solved  using  the 

ju>t 

three  dimensional  Fourier  transforms.  Maxwell’s  Equations  for  the  e  time 
convention  are: 


V  X  E  =  -J“Vi  H 

(i) 

V  <  H  =  j“£  n  +  j 
o 

(2) 

V*H  =  0 

(3) 

<1 

"Ml 

Ml 

II 

(4) 

o 


where  J  represents  the  impressed  electric  current  source  term,  It  will  be 
assumed  that  the  coordinate  axes  have  been  so  oriented  that  €  has  the  form. 


r  ^  '  e  ,  '  A  /  A  a  /  A  A  f  A  A  . 

€  =  xxt  -Jxyc  +  1  y  X  €  +  y  y  €  +  ZZtA 

2 


(5) 


2 


where  x,  'y  and  z  are  unit  vectors  in  the  cartesian  system  and  €q  is  the  free 
space  dielectric  constant. 

It  is  well  known  that  in  a  megneto-ionic  medium  such  a  form  for  €  results 
when  the  z-direction  is  oriented  along  the  DC  magnetic  field. 

As  a  first  step  we  shall  point  out  why  the  conventional  approach  of 
introducing  the  vector  potential  A  through  the  relation 

H  =  V  X  A  (6) 

which  satisfies  (3),  is  not  helpful  in  reducing  the  field  equations  to  a 
simpler  form.  Using  (6)  in  (1)  we  readily  obtain  the  standard  form 

E  =  -j<*V  A  +  V  V  (7) 

and  the  substitution  of  (6)  and  (7)  in  (2)  gives 

V  x  V  X  A  =  V  7>A  -  V2 A  =  w2!!*  f  A  +  j«6  f  V  V  +  j  (8) 

~  o  o 

When  €  is  a  scalar  (O,  (8)  is  readily  simplified  by  introducing  the  relation 

V  V*A  =  j<*>€  C  V  V 
o 

which  reduces  the  equation  for  the  vector  potential  to 

v2!  +  <^u€  e  a  =  -  J 

n  o 

However,  the  corresponding  relation  for  the  anisotropic  case,  viz., 

7  V-A  =  j<*>cjr  v  V  (9) 

is  incompatible  for  a  general  6  as  may  be  easily  verified.  Hence,  as  pointed 
out  earlier,  no  particular  advantage  is  gained  by  introducing  the  vector 
potential  in  the  anisotropic  case,  and  we  shall  find  it  convenient  instead  to 
treat  the  field  equations  directly.  We  now  proceed  to  do  this  in  the  following. 
Elimination  of  H  from  (1)  and  (2)  gives 

VXVXB  =  k2fl  -jWjx  J,  k2  =  £ 


o 


(10) 
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Because  of  the  particular  form  of  F  we  shall  find  it  convenient  to  work 
with  a  different  co-ordinate  system  in  which  the  €  tensor  diagonalizes.  Such 

^  A  A 

a  system  of  co-ordinates  is  defined  below.  Let  the  unit  vectors  U  }  U  ,  U  in 

1  &  <5 

the  new  system,  hereafter  called  the  U-system  be  defined  by 


u1. 


1  ,A  A. 

=  —  (X  +  jy), 

n 


A  A 
(x  -  jy), 


(ll) 


It  may  be  readily  verified  that  the  following  vector  relations  are  satisfied 
A  A  A 

by  Ux,  U2  and  Ug. 


A  A 
U1  •  U1  -  0 


Ui  x  Uf  =  0 


/  A 
U,  .  U  =  1 
1  2 


A  w  A  A 

U1  X  U2  "  ^  U3 


U  .  ft  «  0 
2  2 


A  A 

u  x  u  =  0 
2  2 


61  '  03  *  0 


A  A  A 

°3  *  *1  -  -J  U1 


A 


(12) 
=  1 


The  above  relations  will  be  found  useful  later. 
An  arbitrary  vector  F  given  by 


F 

=  x  F  +9 

F  + 

z  F 

(13) 

X 

y 

z 

can  be  written  in  the 

U-system 

as 

A  A 

A 

(14) 

F  = 

U,  F  +  U 

F  + 

U0  F 

11  2 

2 

3  3 

-  A 

-  A 

A 

where  =  F  .  F^ 

=  F  .  U 

and  =  F 

•  U3 

It  may  be  easily 

varif ied 

that 

Fl  =  jr  <F*  ’  J  V' 


F.  =  ~  (F  +  j  F  ) , 

2  Ji  x  y 


F  =  F 
3  z 


(15) 


Now  the  representation  for  C  in  the  U-system  is 


=  A  A,  A  A  „  A  A  „ 

«  =  U  Un  €,  +  U.  U  eo  +  U  U_  € 

1  21  212  333 


(16) 


where 
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e=€+€/  e  =  e  -  €/  e  =  e 

1  2  3  z 

The  expression  for  the  operator  V  is  as  follows 

V  -  “l  dl  *  \\  *  »3  d3  <17) 

where 

1/8  8\  1/®  8  \  8 

d3=^ 

After  going  through  some  algebra  and  the  use  of  (16)  and  (17),  the  inhomogeneous 
field  equation  for  E  given  in  (10)  may  be  written  in  the  matrix  form  as 


♦  r-<V2  +  V 


d  d 
1  3 


-(dd  +  dd 
1  2  3  2  3 


3  2 


dld3 


"2dld2 


rV 

■v 

E 

J 

i 

2 

2 

- 

Le  - 

L  J  J 

(18) 


The  above  are  the  desired  matrix  equations  for  the  fields.  The  solution  for  (18) 
will  be  obtained  by  taking  the  three-dimensional  Fourier  transform  of  the 
equations  and  by  subsequent  inversion  of  the  transforms.  This  is  discussed  in 
the  next  section. 


3. 


Transform  method  of  Solution 

Let  the  three  dimensional  Fourier  transform  6  of  E  be  defined 
00 


+j(k  it  +  k  y  +  k  z) 

E1  2  3*x>y’z)  e  X  y  Z  dxdydz 


by 


(19) 


and  let  the  transforms  Ck  of  J  ^  _  be  similarly  defined.  Using  the 

^1,2,3  1,2,3 

property  of  the  transforms,  the  differential  operators  d  ,  d  ,  d  may  be 

X  «  J 

written  as 


d 


1,2 


(k  +  Jk  ) 
x  y 


-J 


k 

z 


(20) 
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For  convenience,  we  shall,  introduce  the  polar  form  through  the  definitions. 


k  =  r  sin  vp  cos  a 
x 


k  =  r  sin  +  sin  a 

y 


k  =  F  cos  (21) 
z 


and  rewrite  (20)  as 


sin+ 


+J  « 


d  =  -  J  ^  cos  4> 
3 


(22) 


Taking  the  transform  of  the  matrix  equation  (18)  we  may  then  derive  the 
desired  equation. 


r-k2  . 

r£,  0  0  _ 

.Is  _ 

o 

1 

o  «2  0 

-  0  C  €  - 

L 

V 

3 

r2  .sin  4*  2  |  1  2  |  -2ja 

-*  —  -( — -  +  cos  40  ;  -=  sin  Ye 


i  si„2  +  e2M 


.sin2  Y  2  ,|lV 

■( — -  +  cos  Y)  , 


J~2 


—  sin  4*  cos  Y  e  Ja  ^ 

72 


£ 


—  sin  4*  cos  4*  e 

./2 


j  a 


sin2  4* 


r  6  -j  - j> 


ri 


(23) 


The  determinant  A  of  the  matrix  inside  the  curly  brackets  of  equation  (23), 
obtained  after  working  out  some  algebra  is  expressed  as 


A  =  -  k2  r  (€  sin2  4*  +  €  cos* 
o  3 


40  +  k*  ^  f  €  3  €  ^ 1  ^  cos2  ^  +  (*2  -  **  2)sin2 


k6  (C 2  -  €  2)C 
o  3 


=  -k2  (€  sin2  ^  +  €  cos2 
o  3 


+>  {  r 2 .  „2  <+,]  (r 2  -  „2  <+>', 


(24) 


6 


2  2 

where  n  ,  n  are  the  roots  of  A.  The  roots  are  given  by 
_2  ,2  B(4o  ±  ^(40  -  4A(4>)  c<4>) 

"l,2  *  ko  -  2AW) -  (25) 

where 

A(40  =  €(sin2  4*  4*  €  cos2  40 

w 

B(40  *^€36  (l  +  cos2  41)  +  (€ 2  -  6  2)  sin2  4*j 

c(4»)  =  (€2  -  c/2)  e3 

It  may  be  verified  that  n^  and  are  the  indices  of  refraction  for  the 
extraordinary  and  ordinary  waves,  respectively,  in  the  plane  wave  case.  This 
is  to  be  expected  since  the  condition  that  a  finite  solution  of  (23)  exist 
even  when  A  =  0  is  that^^  =  ^  ^  m  0,  which  obviously  corresponds  to  the 

plane  wave  case.  It  should  also  be  pointed  out  that  the  coefficient  of  r 
in  the  expansion  for  A  is  identically  zero  and  hence  the  third  order  determinant 
in  F2  has  only  two  roots.  The  coefficient  of  F6  ±3  zero  because  of  the 
property  of  the  operator  (V  x  V  x)  which  makes  the  matrix  corresponding  to 
this  operator  a  singular  one. 

It  may  be  verified  that  the  eigenvalue  equation 

|M  -  *  M2|  *  0  (26) 

where  M  and  M  are  (nxn)  matrices,  has  only  (n-1)  finite  roots  if  the  rank 

1  A 

of  M  is  (n-1).  Another  observation  about  A  is  that  it  is  Independent  of  a. 

The  solution  of  the  fields  in  the  transformed  domain  is  now  a  straight¬ 
forward  step  involving  the  Inversion  of  the  matrix  operator  in  (23).  In  the 
matrix  form  the  solution  is  written  as 


1 

=  -J«v 

-au(40/A  «12(4»)/A 

E2 

«21(+)/A  .22(+>/A  «2,<+>/a 

h 

*3. 

_.51<+)/a  .„<+>/a  .33(+)/A. 

i 

3* 

(27) 
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where 


11 


=  —  sin2  4*  -  r  ko  j  C2  Sin2  ^  Sin2  ^  +  cos2  ^ 


+  k  Co  €  (28) 

o  3  2 


■2J OE  ^  2  ij  /p2  2  ,  x 

ai2  - e  r  91"  +  V 


(29) 


a  =  e  —  sin  4*  cos  4*  (T  -  k2  €  ) 

13  7a  02 


(30) 


a2i  =  e2J“  F  811,2  *  <r2  '  k!  V 


(31) 


22  2 


sin2  4<  -  r  k2  (  C  sin2  4*  +  €,(-5  sin2  'l'  +  cos2  40  1  +  k4  C  €  (32) 

°  I  1  32  J  o  1  3 


a  =  e1^  —  sin  4*  cos  4>  (F2  -  k2  C  ) 

23  JZ  01 


(33) 


p4  2 

a  =  e^*  —  sin  4*  cos  4*  (T  -  k2  €  1 
31  /X  o  2' 


(34) 


r2  (T2  _  k2  c  )  sin  4>  cos  + 

32  y2  °  1 


(35) 


«33  =  ^  COS2  4*  -  k2  r2  (6 2  +  (i  Sin2  4*  +  cos2  40  +  k^  6^ 


(36) 


and 


€i=€  +€  > 


«  =  e  -  e‘ 


€  =  €  . 
3  z 


The  main  problem  of  evaluating  the  inverse  transforms  still  remains  and 
we  discuss  this  in  the  following. 


8 


4,  Evaluation  of  Inverse  Transforms 

We  are  now  ready  to  get  into  the  core  of  the  problem,  namely  that  of 
evaluating  the  inverse  transforms  of  the  fields.  It  should  be  pointed  out  that 
the  asymptotic  evaluation  of  the  inverses,  at  distances  far  away  from  the 
sources,  is  a  relatively  simple  task  and  has  been  discussed  by  Bunkin*.  In 
this  paper  we  shall  attempt  to  evaluate  the  integrals  for  arbitrary  distances 
and  our  results  will  be  useful  for  far,  near  and  intermediate  field  calculations. 
To  the  knowledge  of  the  author  a  complete  solution  such  as  this  has  not  been 
discussed  elsewhere. 

Let  us  assume  that  the  current  source  is  an  infinitesimal  one,  i.e. 


J  =  (x  c  +  y"c  +  zC  )6(x  -  x  )  6(y  -  y  )6(Z  -  z  ) 
x  y  z  o  o  o 

=  (VC  +  u  C  +  U  C  )  6(X  .  X  )6(y  -  y  )6(z  -  *  > 
1  1  2  2  3  3  o  o  o 


(37) 


where  ^  is  the  Dirac  &  and  C, .  C  ,  C  and  hence  C  ,  and  C  are  constants, 

1'  y’  z  i9  2  3  9 

and  may  be  readily  obtained  from  the  orientation  and  moment  of  the  dipole. 
For  an  arbitrary  distribution  of  sources  in  space,  an  integration  over  the 
distribution  would  yield  the  desired  result.  For  an  infinitesimal  source  at 
the  origin  the  transform  components  are  readily  obtained  by  the  use  of  the 
property  of  ^  and  are  given  by 


^1,2,3  °1,2,3  (38) 

For  simplicity,  we  shall  carry  out  the  steps  for  evaluating  the  inverse  trans¬ 
forms  for  the  case.  C  =  1  and  C  =  C  =  0.  It  is  obvious  that  through  the 

X  Z  j 

use  of  superposition  we  can  readily  derive  the  complete  solution  when  the 
current  coefficients  C  ,  C  are  non-zero,  by  following  through  similar  routines. 

«  j 

From  (27),  (28)  and  (38)  we  have 


anl<r>  a) 

t  -  n  =  1,2,3  (39) 

V  A(I>) 


The  inverse  transform  of  a  function  ^(r,^, a)  is  given  by 


9 


F  = 


(2ir)' 


-j(k  x  +  k  y  +  k  z) 

(k  ,k  ,k  )e  X  y  Z  dk  dk  dk  (40) 

x  y  z  x  y  z 


-00 


or  in  the  polar  form  as 


oc>  i T  2it 


F  = 


(27 r)‘ 


J J J  (T  *\>  ®  sin  ^  cos  W  +  cos  9  cos  ^ 


0  0  0 


r  sin  4^  da 


(41) 


where  x  =  R  sin  0  cos  <P,  y  =  R  sin  0  sin  <P  and  z  =  R  cos  0. 

Hence  the  inverse  of  (39)  may  be  written  as 

00  IT  2f 

*  ^  -J*!*  I  f  f  J1 i  (T*  ^.“JFR(sin  6  sin  41  cos  (a  -  <P)  +  cos  0  cos 

-*<M>3  ill  a  *  ’a>e 


0  0  0 


•  r  sin  4*  da  dl»  d^ 


(42) 


Let  us  consider  the  case  n  =  1.  From  (28)  we  see  that  a^  is  independent  of 
a  hence  we  can  carry  out  the  integration  in  (42)  with  respect  to  a  giving 

b _ A;  f  f  r*  eJr«  <=°*  e  +  J  (Ife  ,ln  e  ,1„  +>  sin  +  <n>  dr 

<”>  JA  °  (43. 


where  we  have  made  use  of  one  of  the  well  known  integral  representations  for 
the  Bessel  function. 

We  rewrite  (43)  in  the  form 


1  *  *  7*^  y2  1  9> 


k  (2ir) 
o 


(44) 


where 


10 


00  7T 

I (R,  9)  =  k2  J  J  -ii  e-JlRP<0»^)j^(rRq(0  sln  9  ^  dp 
0  0 

p(0,4O  =  cos  0  cos  4^  (45) 

q(0,4O  =  sin  0  sin  4* 

2 

and  A  is  the  Laplacian  operator.  It  was  necessary  to  reorganize  the  integral 
representation  so  as  to  enable  us  to  change  the  order  of  integration,  which  we 
intend  to  do  shortly.  It  should  be  noted  that  since  a^  and  A  are  both  of 
fourth  order  in  T,  a^/A — *  constant  for  large  r  and  the  convergence  of  the 
integral  representation  for  I(R,0)  in  (44)  is  assured  for  all  R  /  0.  We 
anticipate  the  type  of  behavior  I(R,0) — >  0(1/R)  as  R — *0,  so  we  set  out 
to  separate  the  singular  part.  When  this  is  accomplished,  I(R,0)  may  be 
written  as 

I(R,0)  =  I  (R,0)  +  I.  (R, 0)  (46) 

S  7  I  7 

where  I  (R,0)  is  0(1/R)  for  small  R  and  I.  is  finite  for  all  R  and  0.  The 
s  f 

advantage  of  doing  this  is  twofold.  First,  we  shall  find  that  it  is  possible 
to  evaluate  1^  (the  singular  part  of  the  integral)  exactly  and  hence  to 
estimate  the  very  near  field  behavior,  which  is  dominated  by  the  singular 
part.  The  second  reason  is  that  this  procedure  will  make  it  convenient  for 

2 

us  as  we  shall  see  shortly,  to  perform  the  integration  in  the  expression  for  A  If . 
Separation  and  Evaluation  of  the  Singular  Part  I  (R,0): 

Let  us  consider  the  ratio  a^/A.  From  (24)  and  (28)  we  have,  after 
replacing  6  by  €  +  €  and  €  by  €  -  €/ 

X  a 


2  ail 
*o  A 


(1/2)1^  sin2  4*  -  T?  k2  |  (£  -  t/)sin2  4*  +  &3(l/2  sin2  Y  +  cos2  'V)  | 


+  k4  €  (6  -  e') 

o  3 


(Cg  cos2  4*  +  6  Sin2  40  (1^  -  n2)  (1^  -  n2) 


11 


r2r  Sin  4*  2  2V  ,2,*  '  2  ,,, 

1  ~  (n  +  n  )  +  k  (£-€  )  sin  4j 

'-"2  19  o 


♦  ko  8in2  ^  +  COs2  “^€3(6  -if  )- 


2  ail 
ko  — 


1  _ sin2  4> _  T“1  “2  2 _ 

2  (£_  cos2  4^  +  £  sin2  40  (€  cos2  4>  +  €  sin2  4*)  -  n2)  -  n2)" 


2  2  «4n  4* 


(47) 


+n,  n„ 


Now  let  us  define  I  as 
s 


=  -i  7  f  811,8  *  •  8ln^ — 5-  e-Jr»p<e,+>  j  (r,q(e,+)  <ri>  <ir  (48) 

JJ,  «,  co.2  +  ♦  «  ,ln2  + 

U  U  j 


then 


A--I 


7/ 

A  'A 


r2  sin3  4- 


6  cos  4*  +  €  sin  4* 
0  0  3 


rjrRp  J  (T  Rq)  dj>  dr 


11  3  p8 

2  p  %>  ip 


II 

'0  0 


»in  4< 


-JrRp 


£  cos2  4*  +  £  sin2  4*  ° 


J  (Ikq)  <*l>  (49) 


The  integral  in  (49)  permits  an  exact  evaluaticn  as  shown  in  the  appendix. 
Using  (A-4)  we  then  have 


n  i  a  a  c-i/2 

2  '  p^pP^P  ({  ,2  +,  p2j  1/2 

3 


(50) 


It  is  easy  to  verify,  say  by  solving  the  static  problem,  that  V2^  has  the 
correct  order  of  singularity  at  the  origin.  We  also  note  that  the  expression 
for  Ig  is  independent  of  €  ,  the  off-diagonal  term  of  the  6  tensor. 
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As  a  next  step  we  go  on  to  evaluate  V2^  where  If,  the  finite  part  of  I, 
is  given  by  1^  =  I  -  Ig. 

Evaluation  of  V*lf : 

Although  If  itself  is  finite  at  the  origin,  we  shall  find  that  V2^  is 
not  and  it  has  a  singular  part  of  the  type  O(^).  We  shall  now  attempt  to 
separate  this  singularity  from  the  expression  of  V2^  by  the  same  technique 
used  in  the  previous  section.  We  shall  obtain  an  exact  representation  for  the 
singular  part  and  a  reduced  integral  for  the  regular  part  in  V2^ . 

Using  (45),  (47)  and  (48)  we  have  the  representation  for  If. 

f"  r°°  Sfr4oe-jrRp(0>4') 

ffO^e)  =  J  J  y  2  ja  2,  sin*  <51> 

00  (  ni)(  n2) 

where 


r2  r  sin  4*  (  2  +  2y  +  (£  _  €')sin2  ^  +  k2  e  d  sin2  4  +  cos2  4)1 

LZ  1ZO  032  J 


(52) 


S(I>)  = 


.4  *  „  *  \  2  2  sin  41 

-  ko  £3<e  -  e  >  *  “l  "2  —T~ 


C  cos2  4^  +  €  sin2  4* 
o 


Differentiating  under  the  integral  we  obtain,  remembering  p  =  cos  0  cos  4>, 
q  =  sin  0  sin  4* 


r  2  S(T4) 

/r2  2wr  2  2 

(r  -  ni)(F  "  n2 


e-jrRp 


J  (rRq)  sin  4  dr 
o 


(53) 


Using  the  expression  for  S(^40  from 
two  parts  as  follows 
7T  oo 


V2] 


■i  i 


M(4) 


-j  rRp 


6  cos2  4+6  sin2  4 

O 


(52),  break  the  integral  in  (53)  into 


J  (^Rq)  sin  4 
o 


13 


17 


- 2 - ^ ^ - r  e‘jrRPJ  (rRq)  sin  4*  d? 

2  1 1 1 \  .T*  _2.  . ■p 2  2»  O 


o  o  ^ 3  cos  ^  +  6  sin  +)(r  -  n^xr  -  n2> 


(54) 


where 


M(4*)  =  [  -  sin2  4*  (n2  +  n2)/2  +  (€  -  6/)  sin*4  4*  +•  k*4  AO5  sin*4  4>  4-  cos'4  4*)  ] 

12  O  o  0  4 

(55) 


2  *  ,1  _  .  2 


N(^4>)  =  [k4  €  (€  -  €  4  no  sin2  ^  -(n2  +  n2)M(4^)]r2  +  M(4>)  n2  n 

03  212  12  1< 


(56) 


The  first  integral  which  has  the  behavior  of  the  type  0(— )  at  the  origin 
permits  an  exact  evaluation.  The  integration  with  respect  to  F  in  the  second 
integral  may  be  carried  out  using  the  techniques  of  contour  integration,  and  an 
integral  over  a  finite  range  may  be  derived.  The  resulting  finite  range  integral 
is  very  suitable  for  numerical  computations  because  the  Integrand  is  bounded  in 
the  entire  range.  The  detailed  calculation  of  the  first  integral  and  the 
reduction  of  the  second  integral  are  given  in  Appendices  B  and  C,  respectively. 

We  shall  merely  present  the  final  result  which  is 


irk 


A  -  V-  i  •  «e> 


-2irj 


/m,o  -Jn.Rp 

V/  N(n^(40 ,40  e  Jo(n1qR) 

ACF)  7  2  27 

0  "i<ni  -  V 


sin  4< 


(57) 


+  4 


Tt/2  cos  (p/q) 


1TZ2  co 

■nr/2-9  Jo 


N(n1(4'),4*)  sin  ■  nJR(p  -  q  cos  a)J 
M?) 


,  2  2x 

ni(ni  “  n2) 


sin  4>  da  (ft* 


+  P 


(n  ) (R, 6) 
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A(40  =  €  sin2  41  +  €  cos2  4*,  Q(0)  is  given  in  (A-10)  and  P,  N  (R,0) 

'  ^2 

is  obtained  by  interchanging  n  and  n  everywhere  in  the  two  integrals  appearing 
in  (57).  No  further  reduction  of  the  integrals  seems  possible  for  a  general 
choice  of  elements  in  the  €  tensor,  and  numerical  integration  must  be  performed 
to  evaluate  the  regular  part  of  V2^.  Although  we  shall  not  work  out  the 
details,  it  is  nevertheless  of  interest  to  point  out  that  the  two  integrals 
give  rise  to  an  expression  of  the  type 


R  R 


when  the  medium  under  consideration  is  isotropic.  It  is  also  pertinent  to 
mention  here  that  if  one  is  only  interested  in  the  far  fields,  these  can  be 
obtained  by  making  the  asymptotic  approximations  for  large  R  in  the  integrals 
involved. 

The  desired  final  form  for  E^  is  obtained  by  subsituting  (50)  and  (57) 
in  (44).  We  shall  therefore  conclude  our  discussion  of  the  determination  of 
E^  here. 

The  other  components  of  the  electric  field  may  be  evaluated  by  following 
through  a  similar  procedure. 

5.  Comments  and  Discussions 

In  this  paper  we  have  obtained  solutions  of  Maxwell's  Equations  with 
source  terms,  in  a  medium  with  linear  dielectric  properties.  The  very  near 
field  terms  which  have  a  singularity  at  the  origin  have  been  derived  in  an 
exact  form  and  an  integral  representation  for  the  regular  part  has  been  obtained. 
Since  no  restriction  was  put  on  the  nature  of  the  elements  of  the  €  tensor, 
the  solutions  are  valid  for  complex  elements  and  hence  they  are  useful  for 
treating  a  medium  with  finite  losses.  Although  the  details  are  left  out  here, 
it  is  pertinent  to  point  out  that  the  integral  representation  of  the  regular 
part  may  be  asymptotically  evaluated  for  far  field  calculations.  It  is  also 
of  Interest  to  mention  here  that,  where  €  and  €  are  pure  real  and  have 
opposite  signs,  discontinuous  types  of  solutions  result  because  of  the 
modification  of  the  integrals  involved.  This  aspect  of  the  behavior  of 
the  solutions  is  under  current  investigation  at  the  Antenna  Laboratory  by 
Mr.  Keith  Balmain  and  others. 
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APPENDIX  A 


Evaluation  of  the  integral  appearing  in  (49): 
The  integral  to  be  evaluated,  say  L(R,9)  is 


L(R,0) 


-JrR  cos  9  cos  4* 

— - - - - —  J  (^R  sin  9  sin  *1*)  sin  ^ 

C  cos  ^  +  C  sin  ^  ° 


(A-l) 


oo,  irA 

■ll 


^  4^  +  €r2sin2 


63F  C°S 


J  (^R  sin  6  sin  40™sin 
o 


Introducing  a  change  in  variables  through  the  relations. 


we  derive 


y  ss  r  cos  4*, 
z  =  R  cos  0, 

L(E,0)  = 


■JYZ 


U  =  r  sin  4* 
p  =  R  sin  9 

Jo(uP) 

- 2 —  udydu 


(A-2) 


Performing  the  integration  with  respect  tQ  y  we  obtain  through  the  use 
of  contour  integration  methods 


L(R,  9)  =  ir 


4 

(A-3)  is  a  standard  integral  and  is  discussed  in  books  on  Bessel  functions  . 
The  final  form  for  L(R,0)  reads 


j 


«  ^  x 1/2 
~(€/£  )  uz 

e  Jq(i|J>) 


(€  V 


1/2 


du 


(A-3) 
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L(R,e> 


(«Vl/2  [«*3->*2  *  P21 1/2 


(»-•» 


1 _ w _ 

,  1/2  2  a  ,  .2  q.1/2  » 

«  R(t  cos  o  +  €  sin  o) 

u 


Real  («/£3)l/2>  0 


I 

I 

I 

E 


1 


APPENDIX  B 
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Evaluation  of  the  singular  part  of  V2^: 

As  pointed  out  in  connection  with  (54),  the  singular  part  of  V2^  can 
be  evaluated  exactly.  Let  the  singular  part  be  represented  by  K(R,fi)  where 


K(R,0) 


€3  cos2  4^+6  sin2  ^ 


J  (*Rq)  sin  dl* 

o 


(  B—  1 ) 


and  p  =  cos  0  cos  6,  q  =  sin  0  sin  4*. 

As  a  first  step  toward  evaluating  the  integral,  substitute  the  expressions 

for  n  and  n  from  (25)  into  (55)  and  rewrite  M(40  as 
x  £ 


M(40  S  (€  cos  4*  +  €  Sin  40* 

z  u 


X  f  sin4  4M€  -  6  )2  +  sin2  4*  cos2  4*  ^3^3  +  ^ 


/  o 
2€  )  +  26  cos 
J 


(B-2) 
4  A 


6 


Now  multiply  the  denominator  and  the  numerator  of  the  integrand  in  (A-5) 


by  r  and  introduce  the  transformations 

y  s  r  cos  4^  u  =  T  Sin  \\i 


and 


obtaining 


z  s  R  cos  0 


P  =  R  sin  0 


y  mr 

2K(R,0)  =  k2  /  j 

“'O  '-CD 


f  4tc  c>  .2  2  2  ,  ,,  ^2  4} 

|u  (£  -  6  )  +  u  y  J3  £3  +  26  ‘  26  ^  +  ^3  V  J 


2  ,  2.2,  2  2. 

(C  v  +  £u  )  (y  +  u  ) 


x  e"jYZ  J  (ip)  udYdu 


(B-3) 


Evaluate  the  integral  with  respect  to  y  using  the  method  of  contour 
Integration  and  obtain 
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2k‘2K(R,6) 

o 


ir 


e2  +  e'2  +  e2  _  2(cc/  +  e3c  -  e^') 

<€  -  e3)2 


/ 


-uz 

e  Jq  (up)  d 


+  IT’ 


<**  -  *S  +  4€V  -**3*  +££3+£3€/2“36e/2) 

1/2 

2(|~)  C(€-€)2 

3  3 


(B-4) 


-u(C/C3)1/2z 


J  faP  )  du 
o 


C2  +  C/2  -  cc 

_(H)  _ 3 

l2  '  €(€  -  «„> 


7  e-U<€/€3) 


l/2„ 


J  (uP)u  du 
o 


The  first  term  in  (B-4)  is  the  contribution  of  the  pole  at  y  =  -Ju 
whereas  the  second  and  third  terms  are  due  to  the  residue  of  the  double  pole 
at  y  =  -j<6/e3>1/2u. 

Now  rewrite  the  third  term  as 


j 

~  r\ 


-zu(C/63)1/2 

e  Jq(i^))u  du  = 


1/2 

€3  6 
<r')  Xt. 


j 


-zu(C/C3) 


1/2 


J  ( up) du 
o 


(B-5) 


and  evaluate  the  Integrals  involving  the  Bessel  functions  in  the  same  manner 
as  in  Appendix  A,  obtaining  after  some  simplification 


2v~l  k  ~2  K(R 
o 


>e»  “s  { 


e2  +  e'2  +  e2  -  2<ee'  +  e3£  -  c^') 

irrr? 


(B-6) 


€  £2  -  2££0  +  €£2  +  46^'  +  €  2  -  2 

1  3  3  3  3 
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_ 1 _ 

(€  cos2  0  +  C  sin2  0)1^2 

O 


€  e2  +  e  2  -  ee 

\  ^172  (C  -  €3) 


2  a 
COS  0 


.3/2 


(B-6) 


-  i  «<e> 

j  where  Q(0)  is  the  expression  inside  the  curly  brackets. 

i 

1 

I 

1 

1: 


I 
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REDUCTION  OF  REGULAR  PART  OF 


We  shall  discuss  the  reduction  of  the  integral  say  P(R,0),  which  is  the 


regular  part  of  v  1^  given  in  (54).  Let 


P(R,0) 


N(r,40sin  <\>  e~jFRp 
A(4>)  (I*-  n^r2-  n\) 


J  (iRq)  dr  dl* 
o 


(C-l) 


A(40  =  €  sin2  4*  +  €  cos2  4* 

o 


We  note  that  the  highest  power  of  F  in  N(^,40  is  2,  hence  the  integral 
P(R,0)  is  finite  even  when  R 0.  To  evaluate  the  integral,  we  recast  it  in 
the  following  form  after  a  couple  of  simple  changes  of  variables. 


P(R, 6)  =  2 


N(I>)  8in  4*  e~jFRp 

A<*l»)<r®-  njxr2-  n*) 


J  (rRq)dr  dl< 

O 


(C-2) 


The  advantage  of  the  form  of  (C-2)  over  (C-l)  is  that  (C-l)  is  suitable 
for  integration  with  respect  to  ^  using  the  technique  of  contour  integration 
because  of  infinite  limits.  Before  this  Integration  is  performed,  however, 
it  is  convenient  to  replace  the  Bessel  function  by  its  integral  representation 
and  arrive  at  the  form 


P(R,0)  *  | 


/  K  /  2\ 

(A  -  n^)  (A  -  n2) 


2 

+  IT 


HO  '0 


ir/2 

/ 


n(r>).ln  j,  .-jfrtp  -  i  co.  <■> 


a(4o 


»?  <"■  ■!> 


(C-3) 

da  d^  <»1* 


0 
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Figure  1.  Location  of  poles  and  choice  of  contours  for  the  second  integral  in  (C-3) 

Now  interchange  the  order  of  integration  so  that  the  integration  is  first  performed 
with  respect  to  Following  the  usual  procedure,  extend  the  integral  to  a 
infinite  semicircular  contour  over  a  complex  F  plane,  closing  the  contour  in 
the  upper  half  plane  when  the  coefficient  of  ^  in  the  exponent  is  positive  and 
in  the  lower  half  plane  when  otherwise.  Assuming  n  and  n  have  finite  negative 
imaginary  parts,  however  small,  we  see  that  only  the  poles  at  ^  =  n  and  n  are 

X  A 

enclosed  in  the  lower  half  plane  and  their  negative  counterparts  in  the  upper 
half  plane.  Note  that  p  +  q  cos  a  is  always  positive  for  0  <  4^  <  7T/2  and 
0  <  a  <  ir/2  hence  the  first  integral  in  (C-3)  is  always  evaluated  by  closing 
the  contour  in  the  lower  half  plane.  The  evaluation  of  the  integrals  gives 


P(R.  81  -  V  >«>8>  +  V)  <».»> 

X  z 

where  P ,  x  and  P,  x  are  contributions  due  to  the  poles  n„  and  nrt  respectively. 
in^J  in^J  12  ' 

The  expression  for  P,  x  is 


>(V(«,e)  j 


ir/2  ir/2 


N(n1(4'),40sin  +  e_J  1 


-Jn  R(p  +  q  cos  a) 


"ACfT 


,  2  2x 
ni(ni  "  n2> 


da 


(C-4) 


?3 


it/2  ir/2 
•2J  /  / 


(in  4* 


-Jn  R(p  -  q  cos  a) 

K(n  <4»>/b«in  4»  e 


0  cos  (p/q) 


-1, 


A(40 


,  2  2n 

ni(ni  •“  n2} 


d ad\> 


(C-4) 


-  2j 


_*  .  jn  R(p-q  cos  a) 

ir/2  cos  <P/q>  K(ni(4.)^).ln+  e  1 

-4/2-0  '  0 


w 


/  2  2n 

ni(ni  -  V 


da 


and  P  (R  0)  =  P,  . (R,0)  with  n_  and  n  interchanged  everywhere.  The  integrals 
(n  )  ’  (n  )  2  l 

in  (C-4;  may  be  further  reduced.  To  this  end  rearrange  (C-4)  and  derive 


V0*'* 


4j 


+  4 


-jn  Rp 

n2  N(n1(+),4J)  e  1  cos(n1q  cos  a) 

A  (40  n  "  2 

1(  1  "  2* 


sin  4*  da 


7  l 

V 2-e  "o 


cos  (p/q)  N(n^(v|i)4j)  sin  |  n^p  -  q  cos  ojJ 


M40 


nl(n2  -  n2) 


(C-5) 
■sin  4*  da 


Note  that  in  writing  (C-5)  we  have  made  use  of  the  fact  that 


since 


P/q 


F(4%0,a)  da 


0 


p/q 


cos  6  cos  4*  1 

sin  0  sin  4*  ’ 


for  4*  <  ir/2  -  » 


and  the  integral  with  respect  to  a  is  over  a  real  range  only. 

The  first  integral  with  respect  to  a  may  be  reduced  to  a  Bessel  function 
yielding  the  final  representation  for  P^n  ^  as  given  below 
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